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Observations by the Fermi-LAT telescope have uncovered a significant 7 -ray excess toward the 
Milky Way Galactic Center. There has been no detection of a similar signal in the direction of the 
Milky Way dwarf spheroidal galaxies. Additionally, astronomical observations indicate that dwarf 
galaxies and other faint galaxies are less dense than predicted by the simplest cold dark matter 
models. We show that a self-interacting dark matter model with a particle mass of roughly 50 GeV 
annihilating to the mediator responsible for the strong self-interaction can simultaneously explain all 
three observations. The mediator is necessarily unstable and its mass must be below about 100 MeV 
in order to lower densities in faint galaxies. If the mediator decays to electron-positron pairs with 
a cross section on the order of the thermal relic value, then we find that these pairs can up-scatter 
the interstellar radiation field and produce the observed 7 -ray excess. We show that this model is 
compatible with all current constraints and highlight detectable signatures unique to self-interacting 
dark matter models. 


The Galactic Center excess. Recent Fermi-LAT ob¬ 
servations of the Galactic Center (GC) of the Milky Way 
have uncovered a stunning 7 -ray excess compared to 
expectations from diffuse astrophysical emission mn- 
While these studies differ in the astrophysical background 
models, they all agree on three key features of the 7 -ray 
excess: ( 1 ) the spectrum is strongly peaked at an energy 
of approximately 2 GeV, with a low-energy spectrum that 
is harder than expected from 7r°-emission, (2) the excess 
radially extends to at least 10° from the GC, following 
an emission profile that falls with distance (r) from the 
GC as r““ with a = 2.0 - 2.8, and (3) the excess is 
roughly spherically symmetric, without any evidence of 
elongation parallel or perpendicular to the galactic plane. 

While other explanations have been discussed nisi nu¬ 
lls], dark matter remains a compelling possibility. The 
detection of an excess with the same spectrum toward 
dwarf spheroidal galaxies surrounding the Milky Way 
would verify this possibility. However, no equivalent sig¬ 
nal has been observed in dwarf spheroidal galaxies | 16 j . 
which stands in mild tension with some models of the 
GC excess HZ]. 

In fact, dwarf galaxies have long challenged our under¬ 
standing of the nature of dark matter. The dark matter 
halos of dwarf galaxies have constant density cores [181- 
nn, in contrast to the cuspy profile predicted by simula¬ 
tions of cold collisionless dark matter (CDM). Addition¬ 
ally, CDM predicts a population of dwarf halos that are 
systematically denser than the dwarf spheroidal galax¬ 
ies in the Milky Way [52], Andromeda [53], or Local 
Group [211 Hi]. A compelling solution to these chal¬ 
lenges is to assume that dark matter strongly interacts 
with itself |26l [27] . Recent simulations have shown that 
nuclear-scale dark matter self-interaction cross sections 


can produce heat transfer from the hot outer region to 
the cold inner region of dark matter halos, reducing the 
central densities of dwarf galaxies in accordance with ob¬ 
servations [23H3T] . 

Connection to dark matter self-interactions. We 

explore the intriguing possibility that the GC 7 —ray ex¬ 
cess is caused by Inverse Compton (IC) scattering of en¬ 
ergetic e“'"e“ from dark matter annihilation, and the ab¬ 
sence of the GeV 7 -ray signal in dwarf spheroids is a nat¬ 
ural consequence of self-interacting dark matter (SIDM) 
models. Our key observations are as follows: 

• Energetic e“'"e“ from dark matter annihilation (or 
another extended source distribution) can effec¬ 
tively produce 7 -rays in the GC through IC and 
bremsstrahlung, due to the high interstellar radia¬ 
tion field (ISRF) and gas densities in this region. 
The IC emission can explain the peak of the GC 
signal (at 2-3 GeV) for dark matter masses in the 
approximate mass range of 20-60 GeV. The crucial 
requirement is the presence of a new source of e“'"e“ 
with energies larger than 20 GeV, which produce 
7 -rays with peak energy of ^ (20 GeV/me)^i?isRF, 
with typical ISRF photon energy Aisrf ~ 1 eV. 

• The AMS-02 constraint [32] demands a softer elec¬ 
tron spectrum than direct annihilation to e"'"e“ and 
hence annihilation through a light mediator is a 
natural solutiorQ 


^ Another possibility, which we do not explore here, is direct an- 
nihilation to , with e and r channels suppressed. 
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• A nuclear-scale dark matter self-scattering cross 
section requires a dark force carrier with a mass be¬ 
low ~ 100 MeV [55H55] . Annihilations through this 
mediator can kinematically couple only to e“'"e“ 
and neutrinos in the standard model sector. 

• The e“'"e“ produced via dark matter annihilation do 
not produce appreciable 7 -rays from dwarf galaxies 
due to their low starlight and gas densities. The 
dominant signal in these systems should be due to 
final state radiation, with a cross section suppressed 
by ctem- 

• The SIDM density profile in the central region of 
the Milky Way is determined by the bulge poten¬ 
tial [36]. Models of the galactic bulge imply that 
the dark matter density increases to within 1 - 2 ° 
from the GC and the annihilation power is signif¬ 
icantly enhanced compared to the predictions of 
SIDM-only simulations. 

A hidden sector dark matter model. We consider 
a simple hidden sector model in which a 50 GeV dark 
matter particle couples to a vector mediator (j). The 
relic density in this model is set by the annihilation 
process XX 4‘4> with an annihilation cross section 
4.4^ X 10“^® cm^s“^ for a Dirac particle m, where ^ 
is the ratio of the temperature of the hidden sector to 
that of the visible sector at freeze-out |^. We assume 
that the hidden and visible sectors are coupled through 
kinetic [331 or Z-mixing |3D| leading to 0 —>■ e+e“ de¬ 
cays m- The dark matter mass utilized here illustrates 
our main points; the mass could be large or smaller by 
about 50%, depending on the details of electron energy 
loss in the GG. 

In order to compute the secondary emission from this 
model we utilize the software PPPC4DMID |32| , which pro¬ 
vides the solution for one-dimensional diffusion with spa¬ 
tially dependent energy losses. We use the “MED” dif¬ 
fusion parameters listed in PPPC4DMID |32|- This soft¬ 
ware calculates the 7 -ray spectrum from IC scattering as¬ 
suming an interstellar radiation field energy density from 
GALPROP |33|, an exponential magnetic field profile jUj, 
and negligible bremsstrahlung losses, which is a good ap¬ 
proximation for the > 10 GeV electrons under consider¬ 
ation. We tested the PPPC4DMID spectrum by writing an 
independent code that solves the one-dimensional diffu¬ 
sion equation assuming spatially-constant energy losses 
and found good agreement. 

In Fig. [^ we compare the intensity and spectrum 
from our model to the region of interest (ROI) 1 of 
Ref. [TU] , which includes regions within 5° radius (about 
750 pc projected distance at the GG) excluding latitudes 
|5| < 2 °, finding good agreement. 

These results show that the secondary IC emission ef¬ 
fectively reproduces the hard spectral bump at an energy 
of ~2 GeV, and the relatively hard spectrum component 
at energies above 10 GeV observed by Ref. [TU|. The 
hard spectrum component is an important discriminator 



FIG. 1. 7 -ray spectrum from Inverse Compton emission and 
final state radiation produced by annihilation of a 50 GeV 
dark matter particle through a light mediator into e''"e“ fi¬ 
nal state. The spectrum is compared to the Galactic Center 
excess ESj. 


of the dark matter mass, as it absent for masses closer to 
20 GeV. 

We estimate the range of cross sections required to pro¬ 
duce this signal as 0.3 —2 x 10“^®cm^/s, corresponding to 
the SIDM density profiles shown in Fig. and discussed 
next. In order to estimate this cross section range we 
noted (using the density profiles available in PPPC4DMID) 
that the IG signal (shown in Fig. is proportional to 
the J-factor (J = f d£p^(£,fi), where £ = line of sight) 
within 5 degrees of the GC at the 10% accuracy level. 
Therefore, we scale the PPPC4DMID result using the J- 
factors for the SIDM density profiles to obtain the cross 
section range. 

Density profile of SIDM. The cross section estimate 
depends on the dark matter density profile. The expecta¬ 
tion from SIDM-only simulations is that the SIDM den¬ 
sity profile would be essentially constant in this region. 
However, when baryons dominate, as expected in the in¬ 
ner galaxy, it has been shown that the equilibrium SIDM 
density profile tracks the baryonic potential [36]. We 
compute the equilibrium SIDM density profile assuming 
two possibilities for the early (before self-interactions be¬ 
come effective) dark matter density profile following the 
method in Ref. [36]: an NFW profile [45] with scale fac¬ 
tor Tg = 26 kpc [33] and the same profile after adiabatic 
contraction m due to the disk and bulge of the Milky 
Way. The early profiles set the boundary conditions for 
the equilibrium solutions. 

We determine the consistency of this scenario by fit¬ 
ting to the composite galactic rotation curve of Ref. [35] 
shown in the top panel of Fig. [^including a black hole 
of mass 4 x IO^Mq, inner and outer spherical bulges 
with exponential density profiles, an exponential disk 
and a spherical halo. The SIDM halo is computed self- 
consistently from the spherically-averaged stellar distri- 
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FIG. 2. The top panel shows the rotation curve data for the 
Milky Way complied by Ref. |48| and fits described in the 
text. In the panel below, the adopted 7 = 1.2 density profile 
(dashed) is compared to SIDM predictions (shaded), with the 
filled points showing the density profile predicted for SIDM 
without including baryons [29] . 


bution. The general non-spherical equilibrium solution 
m 149) would require non-spherical modeling of the 
bulge, which is beyond the scope of the present work. 

A simple iterative procedure suffices to hnd the best 
fit values after varying the outer bulge and disk parame¬ 
ters. The inner bulge and black hole mass are not varied 
as the SIDM halo profile outside 100 pc is less sensitive 
to them. The SIDM halo has two free parameters: cen¬ 
tral density (po) and central velocity dispersion (cto)- We 
fix CTo = 170 km/s and vary po in the interative proce¬ 
dure described above. Equilibrium solutions also exist 
for slightly different values of cto, but they don’t change 
our key results. 

The SIDM fits and the density proHles are shown in 
Fig.d with the edges of the shaded band arising from 
the two assumptions about the early dark matter profile. 
For comparison, we also show the ht to the composite 
rotation curve with an NFW-like dark matter profile p oc 
r“'^(l + with 7 = 1 . 2 , Tg = 10 kpc and local 

density of 0.3 GeV/cm^. We choose this profile since 
7 = 1.2 is consistent with the GC excess fits [HHSlfTO] and 
it also closely tracks the adibatically contracted profile 
mentioned previously. 

Our inner bulge parameters are consistent with 
Ref. IlH] . The main bulge has an exponential scale radius 
of approximately 0.13 kpc and total mass of 8 —9 x IO^Mq 
depending on the model, again consistent with Ref. [48) . 
The scale radius is significantly smaller than the bulge 
radius found in photometric studies of the bulge [5DH5^ , 
perhaps due to additional bulge structures in the inner 


0.3 kpc, for which there is indirect evidence |53j . For our 
purposes, the present model suffices to convert the ob¬ 
served rotation curve to a baryonic potential, which then 
determines the SIDM density profile. 

The main point of this exercise is to emphasize that 
SIDM predicts high dark matter densities at the Milky 
Way center (and in other baryon-dominated galaxies) 
and to explicitly show that these dark matter densities 
are compatible with current rotation curve data. Thus, 
the SIDM J-factor for annihilation is comparable to or 
larger than the GDM predictions. 

Morphology of the excess. Unlike the case of the 
prompt signal, which is proportional to the annihilation 
power (density squared), the IG morphology is affected 
by diffusion and energy losses. In addition, different 
studies of the GC and inner galaxy isHsi uni have in¬ 
ferred somewhat different spatial templates for the sig¬ 
nal. These facts suggest that signal morphology may be 
an important discriminant but more work is required. 

We note that despite the interaction of the e+e“ flux 
with the non-spherical ISRF, the resultant IC signal 
(using PPPC4DMID) is roughly spherically symmetric (to 
within 10% in the ROI), in agreement with present esti¬ 
mates |5] . Diffusion and energy losses may also moderate 
any asphericity in the annihilation power, suggesting that 
the sphericity of the GeV excess can be accommodated 
in this SIDM scenario. 

In comparison, for a prompt signal to be spherical, one 
must appeal to significant gas cooling in the halo centers 
after all the major mergers jSlj since GDM halos with¬ 
out baryons are highly aspherical in their centers [55] . 
However, adiabatic contraction due to this cooling may 
result in 7 « 2 [56] , which would be inconsistent with the 
GeV excess templates. There seems to be no clear pre¬ 
diction for the morphology of prompt gamma-rays from 
dark matter annihilation. 

AMS-02 constraint. The strongest constraint on 
SIDM models stems from observations of the local 
positron fraction by PAMELA [57] and AMS-02 [58] . 
While leptophilic dark matter was originally posited as 
a solution to the rising positron fraction [5SHST] . these 
models have come into tension with constraints from the 
Fermi-LAT [52], H.E.S.S. [53] and PLANGK [51]. In¬ 
stead, working within the framework that pulsars pro¬ 
duce the rising positron fraction [55ti57] . AMS-02 ob¬ 
servations can set strong limits on the annihilation of 
dark matter to leptophilic final states. For dark matter 
annihilation directly to e“'"e“ pairs, these limits can fall 
below the thermal cross section by two orders of magni¬ 
tude [32]. However, we show below that these constraints 
are relaxed considerably if the e+e” pairs are produced 
in the decay of a mediator because of the softer energy 
spectrum. 

In Fig. we produce Galprop models [55] of the as- 
trophysical electron and positron flux at Earth. We add 
a pulsar component modeled by a power-law injection of 
e+e” pairs with an exponential cutoff [55] . We add a dark 
matter annihilation component assuming the hidden and 
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FIG. 3. The e^e“ spectrum (top) and positron fraction (bot¬ 
tom) for the SIDM model, compared to observations from 
PAMELA and AMS-02, respectively. 


visible sectors are coupled through Z-boson mass mix¬ 
ing, with ^ = 1 and an annihilation cross section to e“'"e“ 
of 4.4^ X 10-2® cm®s-V7 = 6.3 x cm^s"^ [H]. 

This cross section is within the range computed previ¬ 
ously for SIDM models that explain the GC GeV ex¬ 
cess. We use the 7 = 1.2 profile with a local density of 
0.3 GeV cm”® shown in Fig. [ 2 ] since it tracks the up¬ 
per edge of the shaded band near the solar location and 
results in a conservative AMS-02 constraint. Because 
the best-fit Galprop diffusion parameters calculated by 
primary-to-secondary observations |69j may not repre¬ 
sent the average diffusion parameters for leptons, we test 
an ensemble of diffusion parameters and find the model 
producing the best combined fit to the AMS-02 positron 
fraction and the PAMELA e+e” flux. The resulting dif¬ 
fusion parameters are not far from those calculated for 
cosmic-ray nuclei in |69j . Specifically we use a diffusion 
constant of 9.1 x lO^® cm^s”^, a half-scale height of 
6.6 kpc, an Alfven velocity of 30.5 km s”^, and a primary 
cosmic-ray electron spectrum following a broken power- 
law falling as below 11.4 GeV and as at 

higher energies. We adopt charge-dependent solar mod¬ 
ulation, with amplitudes of (f>e+ =171 MV for positrons 
and (/>e- = 54 MV for electrons. 

This result shows that annihilations through a light 
mediator can reproduce the intensity of the GC ex¬ 
cess while remaining consistent with AMS-02 constraints. 
We note that the formal fit for AMS-02 data is quite 
good (x^/d.o.f. = 1.32), but poor for PAMELA data 
(x^/d.o.f. = 3.07). However, we find this to be of a simi¬ 
lar quality to fits produced without a dark matter compo- 
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FIG. 4. Dark matter self-interaction cross section as a fnnc- 
tion of the dark matter relative velocity for mediator mass 
10 MeV (Blue), 50 MeV (red), and 100 MeV (Green). SIDM 
models with ut/ m^ ~ 0.5 — 50 cm^/g on dwarf scales can pro¬ 
duce constant density cores in dwarf galaxies in accordance 
with observations |31 |. 

nent. Interestingly, updated measurements by AMS-02 
offer the exciting possibility of constraining or detecting 
SIDM annihilation. We note that these findings are con¬ 
sistent with [32], given that we adopt a local density of 
0.3 GeV cm”®, as opposed to 0.4 GeV cm”®, and noting 
that the annihilations through a light mediator soften the 
resulting e"'"e- injection spectrum, making it more com¬ 
parable to annihilation through than to direction 

annihilation to e"'"e- pairs. 

SIDM solution to small-scale structure formation 
problems. In Fig. [^ we show the dark matter self¬ 
interaction cross section as a function of the dark matter 
relative velocity. For 50 GeV dark matter, the mediator 
must be less than 50 MeV for self-interactions to solve 
anomalies on dwarf scales. It is interesting to note that 
for attractive interactions ct is enhanced when the medi¬ 
ator mass is 50 MeV because of the s-wave resonance [70] . 
The cross section ctt drops slightly on Milky Way scales, 
but is still large enough to effect the Milky Way halo [35] ■ 
On cluster scales, dark matter self-scattering is highly 
suppressed as l/r;^, because the momentum transfer is 
much larger than the mediator mass and dark matter 
self-scattering occurs in the Rutherford scattering limit. 
Therefore, the model is fully consistent with constraints 
otI^x ^ 1 cm^/g for v ^ 3000 km/s from the Bullet 
Cluster m and cluster shape constraints m 
Other detectable features. In addition to producing 
a correlation between the density profile of dark matter 
and the galactic bulge, as well as producing a signifi¬ 
cant contribution to the AMS-02 positron fraction, the 
dark matter model described here may be tested through 
radio observations. Specifically, the large e^e“ flux pre¬ 
dicted by our model may be able to explain (or be con¬ 
strained by) the Green Bank Telescope radio continuum 
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observations towards the GC m, the isotropic emission 
detected by ARCADE-2 [74U77] and the observation of 
hard-spectrum radio filaments in the GC m- 

Additionally, the proposed SIDM particle can scat¬ 
ter with nuclei through the mediator (j), leading to di¬ 
rect detection signals [41]. Since the mediator mass is 
comparable to the momentum transfer of nuclear recoils, 
the event spectrum has a non-trivial momentum depen¬ 
dence miiig, which could provide a smoking-gun signa¬ 
ture for SIDM. 

Summary. We show that the GC excess can be ex¬ 
plained through secondary emission from e“'"e“ pairs pro¬ 
duced in dark matter annihilation events, a scenario 
which naturally predicts suppressed 7 -ray emission from 
dwarf spheroidal galaxies. This class of models is well- 
motivated in the context of SIDM models posited to ex¬ 
plain anomalies in the dark matter density profiles of 


dwarf spheroidal galaxies and other baryon-poor galax¬ 
ies. These models make unique predictions, which could 
be tested in the near future. 

While this paper was in preparation, several related 
papers were submitted by other groups [80111I| ■ We note 
that our favored SIDM mass range is consistent with the 
analysis of light-mediator models by [50] ■ 
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